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WAVE  GENERATION  BY  POINT  SOURCES 

by 

W.  A.  Nierenberg 


I.  Introduction 


The  problem  of  wave  generation  In  water  of  arbitrary  depth 
by  complicated  situations  such  as  moving  ships,  surface  disturb- 
ances, etc.,  has  been  dlsousned  many  times.  However,  the  treat- 
ment of  the  simple  model  of  a point  souroe  at  a fixed  olroular 
frequency  has  not  been  published.  The  advantages  of  suoh  a model 
are  obvious.  The  judicious  disposition  of  suoh  souroe a in  space 
can  be  used  as  a rough  model  for  almost  any  souroe  of  disturbance. 
The  Fourier  synthesis  of  the  solution  oan  be  used  to  interpret  the 
disturbance  by  a souroe  varying  with  the  time  in  any  manner. 

In  this  report,  the  exact-  solution  for  all  distances  from  the 
souroe  for  a point  source  in  a linearized  theory  is  obtained  a la 
Lamb.  Application  is  made  to  a dipole  source  whose  dimension  is 
short  compared  to  the  wavelength  of  the  generated  disturbance  at 
large  distances. 


II.  The  Differential  Equation  and  the  Boundary  Conditions 


In  this  problem,  the  equation  to  be  solved  is  the  equation 
of  continuity  on  the  velocity  potential  <£  . 

V2  ^ = O (i, 


The  solution  to  this  equation  will  be  made  to  satisfy  the  free 
surface  boundary  condition  and  the  condition  of  a spherical  source 
at  a depth  h from  the  surface.  The  It  -axis  coincides  with  the 
vertical,  with  ^ being  measured  as  positive  in  the  upward  direc- 
tion. The  water  depth  is  H.  The  usual  technique  of  splitting 
the  region  into  two  parts  using  a horizontal  plane  through  It-  3(i 
will  be  followed.  The  region  o will  be  labeled 

with  the  subscript  I,  the  region  with  the  sub- 

script II.  The  elementary  solutions  for  tne  two  regions  are: 
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§x  = eiwt(Aek*+Be-k*)Jo0 »/«.) 
# = e‘at(C  +De'^)J,(kn) 


CO  la  the  circular  frequency  of  the  harmonic  source,  k Is 
the  separation  constant,  A,  B,  C,  D are  arbitrary  constants,  and 
r is  the  cylindrical  radius  from  the  origin.  The  boundary  con- 
ditions to  be  satisfied  are: 


^3^  ^ at  ’ % 2/ 


2-t  w 


3 = 3i 


'c>jjt>Tr 


Equation  4 is  the  well  known  free  surface  condition.  In  Eq.  5 

refers  to  the  "ft  -factor  of  Eqs*  2 and  3.  The  disconti- 
nuity in  the  component  of  the  velocity  is  illusory,  for 
the  final  linear  combination  will  involve  an  Integral  like 


k- r 
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f <TcC kit)  kith  whioh  is  equivalent  to  <Tc/ X)  . Equation  6 

represents  continuity  in  pressure.  Equation  7 requires  vanishing 
vertical  motion  at  the  bottom  boundary.  V is  a constant  whioh 
will  be  interpreted  subsequently. 

The  indicated  substitutions  are  performed.  The  constants  are 
evaluated.  Finally,  the  simple  linear  combination  of  summation 
over  all  k is  chosen,  yielding 


f = ^ e:atf . — 


J J0Ck/i)e£j^ 


There  is  still  an  ambiguity  in  the  path  of  integration:  the  inte- 

grand has  two  pole 8 on  the  real  axis  given  by 


/>o tuL  k H 


(10) 


These  roots  are  of  the  same  absolute  value  but  opposite  in  sign 
We  will  oail  them  I k0.  The  path  of  integration  is  chosen  to 
avoid  the  pole  by  a small  semicircular  excursion  into  the  upper 
half  of  the  complex  plane*  (Fig*  1)  This  will  guarantee 


Figure  1 
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that  the  final  solution  will  consist  only  of  outgoing  waves  at 
infinity.  We  note  that  the  real  solutions  of  Eq.  10  give  the  cor- 
reot  phase  velooity  for  the  outgoing  waves. 

Finally,  there  ia  the  question  of  verifying  the  singularity 
at  the  origin  and  obtaining  the  normalization  of  the  source.  In- 
spection of  Eqs.  8 and  9 shows  that  the  Integrals  go  to  infinity 
as  both  and  « Furthermore,  the  infinite 

contributions  come  from  the  very  large  values  of  k in  the  inte- 
grand. Therefore,  to  th8  approximation  that  ft-5>o  and 

we  may  approximate  the  solution  by 


(ID 


& e‘wt/e  hL^Z ;cto)4k  = 4 K 


(12) 


R = 


It  follows  that,  in  the  neighborhood  of  the  singularity,  the 
solution  is  represented  by  a source  whose  peak  strength  is  V 
units  of  volume  per  second. 


III.  The  Evaluation  of  the  Integral 


In  addition  to  the  two  real  poles,  the  integrands  of  Eqs.  8 
and  9 have  an  infinite  number  of  pure  imaginary  poles,  symmetrical 
as  to  sign,  given  by  Eq.  10.  These  roots  will  be  denoted  by  kn, 

/K/  * il,  18,  13.  There  are  two  limiting  oases  of  general  in- 
terest. 
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Case  1. 


or  the  shallow  water  approximation. 


In  the  limit  we  finds 


/ Hk„  I - ^ 


1,2,3 (13) 


Case  2. 


There  is  no  simple  approximation 


but  the  early  roots  start  as 


n ■ 1,2,3. 


(14) 


For  large  values  of  n,  the  solution  turns  over  to  agree  with  (13) 
but  the  first  roots  are  the  only  practical  ones  in  the  limiting 
solution.  If  more  than  a few  roots  are  required,  the  series  ex- 
pansion from  the  origin  is  more  useful.  Since  there  are  no  branch 
points,  the  integral  is  evaluated  by  elementary  technique  yielding 


_ - ,tuit  -«> 

> 1 ” 


I +•  X hi 


(15) 


where  H§  ( U,  ) is  the  zero-order  Hankel  function  of  the  second 
kind.  This  is  the  solution  to  the  problem.  It  is  to  be  noticed 
that  the  artificial  division  of  the  solution  valid  for  two  regions 
has  disappeared  in  the  final  expression  in  Eq.  15.  This  must  be 
so  and  the  fact  that  this  occurs  is  a oheok  on  the  solution.  It 
will  be  convenient  to  regard  Eq.  15  as  the  Green's  l’unctlon  for  the 
problem. 


In  any  practical  application,  this  solution  can  be  used  to 
Idealize  the  situation  by  appropriate  combinations  of  sources  and 
sinks.  As  an  Important  example,  we  shall  consider,  in  addition 
to  the  source  at  the  origin,  another  source  whose  coordinate  is 


S'c)^  i ^ ) ‘-f  “ ° which  is  oscillating  just  7 r 

radians  out  of  phase  with  the  first.  In  this  case,  we  have  the 
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difference  of  the  two  solutions.  Since  the  various  addition  theo- 
rems are  to  he  applied,  the  solution  is  manipulated  in  a prelimin- 
ary way: 


r*~-  o 
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Veiut  3 ,ct*Lk,»(H+fj &nA-L(H+*,-J)  r, 
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(16) 

oWj 


where  is  the  oylindrioal  coordinate  measured  from  the  seo- 

ond  source. 

The  two  terms  of  Eq.  16  may  be  written  in  a more  familiar  form 
through  the  use  of  addition  theorems. 
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This  solution  is  exact.  In  the  next  section  various  approxima- 
tions are  considered. 
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The  Various 


sroximat  lone 


The  Asymptotic 


iroxlmatlon 


In  this  approximation,  the  concern  Is  with  the  wave  am- 
plitude at  distances  such  that  /?« /L  » 1 . The 

customary  technique  Is  to  employ  the  asymptotic  expansions 
of  the  Bessel  functions.  In  addition,  we  recognize  ♦hat  the 
terms  for  n » 0 yield  the  outgoing  waves,  whereae  the  terms 
n * -1.  -2,  eto.,  yield  exponentially  damped  waves  with  dis- 
tance. In  praotlce,  the  damped  terms  may  generally  be  neg- 
lected at  about  one  wavelength  from  the  source.  This  approxi- 
mation gives 

^ oo  • (18) 


In  this  approximation,  the  departure  is  Eq.  18  with  the 
additional  assumption  that  Qo  U<<  1 • The  ex~‘ 

presslon  further  simplifies  to 

jr.  V e:ut-i*'n  -Jf_  v 

/ fwi?,  a 
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This  approximation  Is  Intentionally  Inconsistent  In  the  terms 
of  the  second  order.  The  effect  of  the  vertical  separation 
of  sources  vanishes  In  first  order  and  the  second  order  terms 
have  been  kept  for  the  special  case  d * 0.  It  Is  also  useful 
to  have  an  expression  for  an  isolated  source  to  the  same  ap- 
proximation. Reference  to  Eq.  16  yields 


5>  = ~ r \J 

i n < j nr<P/ l. 


(20) 


3.  The  Dipole  Approximation  in  Shallow  Water 

In  this  approximation,  the  assumption  is  made  that  the 
two  sources  are  considerably  less  than  a wavelength  apazt  or 

that  both  asiJL  /A> 

To  terms  of  the  first  order  only. 
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(21) 


In  effect,  the  single  source  is  reduced  by  the  fraction 
A©  f*/  j f . The  effect  of  the  vertical  separation 

I “ ^ v-  in  shallow  water  is  of  higher  order  and  is 
evaluated  in  Eq.  19. 

There  are  further  useful  approximations  but  they  will 
not  be  listed  here  since  they  are  of  secondary  importance. 


V 


Numerical  Estimate b 


As  an  example  of 
we  consider  the  case 


Since 


the  order  of  magnitude  of  the  effects  Involved, 
/y=  l(K>\  a)  - IT  ho  /ULc  , 9lnele  source. 


t,  ~ K_ 


(20a) 


This  corresponds  to  a wavelength  of  1120  ft.  If  an  overpressure 
of  1 Inch  of  water  (approx.  6 lbs/ft2)  is  considered  at  a distance 
of  1000  ft,  the  peak  flow  volume  In  ft3  per  sec  Is 

V - 9.0  x 103  ft3/sec  or 

V ■ 5.4  x 105  c.f.m. 


If  the  source  Is  a horizontal  doublet  of  spacing  20  ft 
the  V required  would  be  6 x 10®  fts/min  which  Is  considerable. 
For  purposes  of  computation,  It  Is  convenient  to  write  Ea.  20a 
In  terms  of  oO  and  //  alone. 
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